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УДК 51: 621. 891 + 06 ООТ 10.12737/6808 
Модель стратифицированной смазки упругодеформированного радиального подшипника* 
С. В. Митрофанов, Б. Е. Копотун 


На основе уравнений Навье-Стокса для случая «тонкого слоя», уравнения неразрывности, уравнений Ламе и урав- 
нения Дарси, с учетом граничных условий на поверхности шипа и подшипника, на границе раздела слоев, приво- 
дится метод формирования точного автомодельного решения задачи гидродинамического расчета упругодефор- 
мированного радиального подшипника, работающего на двухслойной смазке и обладающего демпфирующими 
свойствами. В результате решения задачи найдено поле скоростей и давлений в смазочных слоях и получены 
аналитические выражения для основных рабочих характеристик упругодеформируемого радиального подшипника. 
Проведен численный анализ полученных аналитических выражений для безразмерной несущей способности упру- 
годеформируемого радиального подшипника с пористым слоем на поверхности шипа. Приведена оценка влияния 
значений параметров, обуславливающих подшипнику повышенную несущую способность. 

Ключевые слова: улругодеформированный радиальный подшипник, стратифицированное течение, адаптирован- 
ный профиль, пористый слой, упругогидродинамический параметр, несущая способность. 


Введение. Имеется практический и теоретический интерес к исследованию упругодеформируемых 
подшипников [1-3]. Выполненные исследования в этом направлении [4, 5] показывают, что упругоде- 
формируемые подшипники обеспечивают большую устойчивость при работе на высоких скоростях и 
малых нагрузках, чем жесткие подшипники. Существенный недостаток существующих расчетных моде- 
лей упругодеформированных подшипников состоит в том, что в них не учитываются особенности взаи- 
модействия смазочной жидкости с твердой круговой опорной поверхностью подшипника, в результате 
которого происходит расслоение смазки на слои с разной вязкостью. В известных работах [6, 7], посвя- 
щенных расчету подшипников, работающих на двухслойных смазочных материалах, опорная поверх- 
ность предполагается жесткой и, кроме того, поверхность подшипника считается сплошной. В работе 
[8] приведена расчетная модель упорного подшипника скольжения, работающего на трехслойной сма- 
зочной композиции с учетом деформации его опорной поверхности, где рассматриваемый подшипник 
не обладает демпфирующими свойствами. Целью данной работы является обобщение предложенного 
в работе [8] метода для случая радиального подшипника, работающего на двухслойной смазочной ком- 
позиции и обладающего повышенной несущей способностью и демпфирующими свойствами. 

Постановка задачи. Рассматривается установившееся стратифицированное течение двухслойной 
смазки в зазоре упругодеформируемого радиального подшипника скольжения, обладающего демпфи- 
рующими свойствами. Предполагается, что подшипник с адаптированным профилем опорной поверх- 
ности содержит упругий слой и неподвижен, а шип с пористым слоем на его рабочей поверхности вра- 
щается с угловой скоростью ©. В полярной системе координат уравнения контуров 9 о ие 


(рис. 1) можно записать в виде: 
СГ =Ю; бо: =к+И; С, :Г'= 0 +Н+6а + ае с0$60 — аАзт 0 + ^'Ё(6); 


С. :/'=Б +е6050 -А$тоб; С, : /'=^ +ес0$0- Азп 0 + ^'/ (0); (1) 


С, :Г'’= г, +ес0$0-— Азпт 0, 
где ое [0,1]; М — толщина пористого слоя, мм; г, — радиус шипа, мм; в=/ -/; Аи о — соответ- 
ственно, амплитуда и частота контурного возмущения в мм и с! (характеризует степень отклонения 





* Работа выполнена в рамках инициативной НИР. 
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контура опорной поверхности подшипника от кругового); ^’/(0) — ограниченная функция, харак- 
теризующая деформацию упругого слоя на поверхности подшипника, мм; г. - г, — толщина упру- 
гого слоя, где г, — радиус внутреннего контура подшипника, прилегающего к смазочному слою, 


мм; г, — радиус деформированного контура упругого слоя подшипника, мм. 


у! _ С 




















х| 




















Рис. 1. Схематическое изображение шипа с пористым слоем в упругодеформированном радиальном подшипнике: 
С, — внутренний контур пористого слоя шипа, прилегающий к непроницаемой поверхности; С, — внешний контур пористого 


слоя шипа, прилегающий к смазочному слою; С, — граница раздела 2-х смазочных слоев; 
С, — внутренний контур подшипника, прилегающий к смазочному слою; 


С, — деформируемый контур упругого слоя подшипника; С, — внешний контур подшипник 


Решение задачи. В качестве основных уравнений примем безразмерную систему уравнений движения 
вязкой несжимаемой жидкости для случая «тонкого слоя», уравнение неразрывности, уравнение Дарси 
и уравнения Ламе: 














0’, ар ди д, 0Р 10Р ТОР 
В 
ог 40 ог 09 05° С0 С’ 00 
0?и, д?и 
—=0, = =0. (2) 
ог ог 
Величины /’, и’, о’, р’ в смазочном слое связаны с безразмерными г, ‚о, и р, соотношениями: 
' , , , ы * шО(и. +Н)}? 
м=к+Нчбг, 9, =Олу,, и! =Оби,, р'=р,р,, р; = 52 НР 
где ш,, о, — компоненты вектора скорости, мм/с; р, — гидродинамическое давление в смазочных 


слоях, Па; и, — динамический коэффициент вязкости, Па:с. 
В пористом слое переход к безразмерным переменным осуществляется по формулам 
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О (3) 
где Р’ — гидродинамическое давление в пористом слое. 
В упругом слое переход к безразмерным переменным осуществляется по формулам 
= +9Г”, ЕЛ 


, * ’ * 
ТЕДА, М, = ММ, Ш =ИЩ,, (4) 


где и’, и — компоненты вектора перемещений. 
Граничные условия на поверхности шипа и подшипника записываются в виде: 














ОР .. ОР 
и =-\№М—- , 20 — 1, 0) = 2п) = 5! = 0; =Р , 
,| г=0 96 ыы о, и р‚( ) р‚( п) РВ. 06 ы р, | 
и, --вв)= 0, 55 | „-ве) = 0, Рь(0) = р, (2) = Р,. (5) 
На границе раздела слоев граничные условия записываются в виде: 
0; до и 
и = ' гай: = о Я ай о ей, р’(0 р 
1 | гов 2 | о, ре о |рие- др И-ай — Гор |ичай о, ел ай’(0) 
| КК +Н)? 7 р 
В(0) =1+ 1050 —т, $п об + ^./ (0), п=е/5, п, =А/5, М= и ) ИВ, =>. (6) 


Граничные условия (5) означают прилипание смазки к поверхности шипа и подшипника, пери- 
одичность гидродинамического давления, а также то, что при переходе через пористую поверхность 
давление меняется непрерывно, а нормальная составляющая скорости определяется законом Дарси. 
Кроме того, на непроницаемой поверхности шипа нормальная составляющая скорости равна нулю. 

Условия (6) означают равенство скоростей, касательных и нормальных напряжений на границе 
раздела слоев, а также условие существования слоистого течения смазки, т.е. требуется, чтобы ско- 
рость точек границы раздела слоев в каждой точке была направлена по касательной к контуру раздела 
слоев. 

Граничные условия для уравнений Ламе записываются в виде: 


ди С 
г -Р, и, 0, (7) 


— | * = и. 
дг” г =^ (0) 


и-ь (0) 








и-ь (0) = 
где р=тахр,, 0 =[0,2*]; р, — безразмерное гидродинамическое давление, найденное в работе [9], 
в случае подшипника с жесткой опорной поверхностью, работающего на двухслойной смазке; 

_ С +а’)и"5* , 
(1-а”)иье(л, + #25, ' 


гидродинамический параметр; 





а’ — постоянная Мусхелишвили; С. — модуль сдвига; ЛМ — упруго- 


В, (0) = 5 0050- т 6, 1,(0)=1+ = 050 - ат соб. (8) 


1 1 1 1 
Полагая толщину пористого слоя достаточно малой, уравнение Дарси усредним по толщине 
этого слоя [10] 


7+1 2 2 
Е а5=0. (9) 
06° 606 С 09 
Точное автомодельное решение системы уравнений (2), удовлетворяющее граничным условиям 
(5) и (6), с учетом (9), ищется в виде: 
Е ду; 
00 


у 





+Ие,®), о, = +И(г,0), 9, =9 9, 
ооо МЕНЕ в 
' ' ' ' в’ 0 № 12’ 00 № № 
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Р=А(6) (5-9) (6-(+1)) +с" (6-7) В(6С- (1+1)1+ р, 


В -.Б 
Е" +). (10) 


Подставляя (10) в (2) и в граничные условия (5) и (6), с учетом (9), получим: 


от я 


у" =&, 5" =&, +86, =0, $," =6, 5, =&, +8 =0; 
$, (0) = 0, 1 (0) = -М№", 5, (0) =1, > (1) = 0, 0 (1) =0, (11) 
5. (1) = 0, у! (©) — > (о), 6, (©) = 6. (а), 0, (©) = 0.(а), 


5” И? -, =! и " [бо 
5: (а) = —2 6, (а), $: (а) = —лу5(а), р, = 0 
Ш и: 1 


[е 1 
[5 (®) + [5,(©)4Е = -№" =В". (12) 
0 а 

Учитывая, что расслоение смазки на слои происходит вблизи неподвижной твердой поверхно- 
сти, т.е. при значениях оа,, близких к единице, условие раздельного течения смазки 
(2, (а) / 5,(о)) = а’(0) в принятом нами приближении удовлетворяются. Из граничного условия 














й, (а) + аб, (а) + | 5,(В)4Е=0 (13) 
следует, что . 
со о СЕ : (14) 
Используя теорему о среднем значении, получим: 
сы к Ре а о с. (15) 
Так как 5, (“=”) < 6,(<), (1-я) <<1, следовательно, с точностью до членов [7 (1- В] 
будем иметь | 
[5 (94 + М" =0, [50 = 0. (16) 


0 а 
Решение задачи (11)-(12) находится непосредственным интегрированием. В результате 


Я 5 Е 


+СЕ+С,, => +СЕ+С,, 





—: 2 
В 2 Е = Е 
о 9 О РО 
и _ Е 2 й Е? Е? 
а а. + Со, о ы 
р =&.,(0) + &7.(0)+с,, р, =&7(0)+ 61 (0)+с., 
0 
90 
2. (0) = | =. (17) 
о (1+1С0$0 —\, $ ®0) + ^./(0) 


Для определения постоянных с, (/ = 2,3.,...13), 6,, 6», си, ©, решим следующую алгебраическую 
систему из 16-ти уравнений с 16-ю неизвестными: 


169 


Вестник ДГТУ. 2014. Т. 14, №4 (79) 
































©х=1 с,=В», 6:=0, С,=Б, С =Р,, 
& 3-6 5+Си 0, & 5+6 +6, =0, 5+4 + С 50.96, т! 
= ы - с (2^) = | 2 | 
© Ро а в) 22 ИЕ. че О С,4+С ей 
а? О а? = а? 
с. + С: + С, г. с. -— С; =0, а +а+с, 4-6 6, =0, 
5 о а? = 90° ? = а 
а ее #5 О 2 +с, =-№ =В. (18) 
Решение системы (18) сводится к решению следующего матричного уравнения: 
М.Х=Ь, (19) 
где х = {С;сисусь:с,}; 6 = {0;0; ба — 68";0;—2}; 
дат 2 2 0 0 
7, (2п) 
1 0 0 2 2 
М=| Ка +1 0 0 ЗКа? -3а? +3 б-ба|. 
5, (2 
Л 0 0 
7, (2п) 
о’(К-1) 0 0 20(К -1) —2 





Решая систему (18), получим 





6+2В° + бКа? — ба? — 2аВ” + аКВ” 
А 7’ 





ё= 
_ 4, (2) 1 
“ 0) (ок-а+1) 
+3К?а* + а АВ’ — 2а?КВ” +В’Ка? — 6Ка^), 
4) , 1 
ь 7, (2п) (ак-а+1) 
+6Ко? — ЗК + Заз К? — ВК -—2В"Ка? + За КВ” — За? К? — ВК? — 6Ка? + За К), 
_4- В’? — 4а? +В’Ка? + 4Ко? +В” 
ИА ’ 
р 4а? + 4Ко? — ЗКо? + В’Ка” + За? — В’? + аВ" — аАВ” +1 
9 А 7 
с, =КС., с, = КС., 





и (3—6? —В'м? + За +В’? — ав” + о А?В” + 








т а(—3За.? — За — 2 В" + За? +В’а? +В +В Ка? +3+ 








С 














А =-—4а° +1+а* — бКо” +4Ко? + Ка“ +4Ка — 2Ка" — 4 + ба”, 
Е (20) 
7, (2п) 


Таким образом, для определения поля скоростей и давлений необходимо найти функцию ^./(6) 





С, 


‚ характеризующую деформацию опорной поверхности подшипника. 


Для определения этой функции воспользуемся приближенной формулой 
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- (21) 


|. (6) И не г 





Таким образом, из найденного выражения для функции ^./(0) следует, что при М —> о, и, 0 


(это соответствует случаю жесткой опорной поверхности подшипника). 
Воспользовавшись формулой (21) выражение для /(0), можно преобразовать к виду 


В(е) = (1+2.^(0)) (1+1 с0$0 — Я, $по0) = 





- [1+ + лсово- А,ао), м = Нее й, = Мы. (22) 
М р 
1+ — 1+ — 


< 


С учетом (22) на основании граничного условия р, (0) = р,(2*) для с,, получим следующее вы- 


ражение: 








с. ни | р, (соваль 1 (23) 
Вы ыы 
М 


Используя (23) для гидродинамического давления Р ‚ окончательно будем иметь 


ры О. ООО. ООО Ра, (24) 


О: м2 >72 р 
Вы ое. А 
М М М 


Основные рабочие характеристики подшипника запишутся следующим образом: 


К, я й = 3 к 
- Бра: [ий й, |5 1721 -1 , 605(@ + 1)2п :| 


Як к. 2% ®-—1 ®+1 























РП. 00 


К, _ а РИ ть ок, 
М) 2° Ф-—1 Ф-+1 


=! 


г Р [6] 2т у" ее 2 Ц 
[ = ты а 1 = _ Ч г . 
) —=] |. т | С, | го (60$ 2л® — 1] +с, | г: (с0$ 2ло 9 





но 
ь С 1 
г. 5 Е 
[ен 


Анализ полученного аналитического выражения для безразмерной несущей способности КА, 


(25) 











(с05 2 5, С 1+ й 


о) 5 (с05 2ло — о 


м 


приведен рис. 2—4. 
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Рис. 3. Зависимость безразмерной 


несущей способности от параметров 
хим. 





Рис. 2. Зависимость безразмерной 
несущей способности от параметров 

о и = . © 
г 1) Р. =0, 01, 
1 М=100, В=-0,1; м 
2) М=1000, В’ =-0,1; 2) я * =—0,1; 


3) М>«®, В=0. 


К у 
г — 


2% 





100 М 


80 


60 





20 


Рис. 4 Зависимость безразмерной несущей способности от параметра Л/ 








ни: Е . 
1 0, 0; 2 0,02, =—0,1; 3) — =0,01, =-0,1. 
ри -0, В и в и в 
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Выводы. Результаты численного анализа, приведенные на рис. 2—4 показывают: 

1. Несущая способность подшипника существенно зависит от значения упругогидродинамиче- 
ского параметра /Л/. На рис. 4 проиллюстрировано, как с увеличением значения параметра Л7 повы- 
шается несущая способность подшипника. 

2. При М > ® несущая способность подшипника стремится к соответствующему значению для 
подшипника с жесткой опорной поверхностью. 

3. С увеличением значения параметра В" Г обусловленного наличием пористого слоя на поверх- 


ности шипа, несущая способность подшипника незначительно снижается. 
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МОБЕЕ ОЕ $ТВАТТЕТЕО ГОВЕТСАТТОМ ОЕ ЕТАЗТТСАН-У ОЕРОВМЕО ВАОТА! ВЕАВ1МС* 


$. У. МИгоГапом, В. Е. Коровип 


Оп пе Баз/5 о те Маиег-51юокез едиайоп$ Юг {те сазе ога 7 /вуег, те сопёпийу едиавоп, [ате едиавоп5, апа Рагсу 
едиавоп, сопздейпоа ве Боипаагу сопайоп5 оп пе рп апа Беагпд !асе аЁ {пе (вуег (тепасе, а тетод оЁ депегайпд ап 
ехасЁ зе йтИаг 5оивоп ю (пе рго ет о! сасшата йе пуагодупати/с евзЁсайу деюгтеад га! Беата иогКпа оп а 
вло-уег ибтсавоп апа ваутпд аатрд ргорегиез 15 ргоиаЕеад. Тре ргоет зо/ивоп (5 гези!ед п епёу/па пе иеосйу 
апа ргеззиге пе4 п пе ибисайпад (вуег, апа п обатпа {йе апа/уйса! ехргезоп5 г те Бах регоптапсе о! те ей- 
Бса/у аеогтеа га@а/ Беагпд. А питепса! апа/уз5 оЁ те оатеа апа/уйса! ехргеззюпз Гог те поп-а@тепзюпа! /сас- 
сатгупад сарасйу о! {те е/азйсайу деогтеа га! Беата и а рогоиз (ауег оп пе рп зигТасе 15 ретогтеа, Тре еТесЕ о! 
Ре рагатег иаиез рота пе Веайпд ий ап псгеазед /саа-сатутпод сарас!!у (5 а5ес5ед, 

Кеуигогаб: е/а°Ёса/у дЕогтесд гад! Беагпа, гаййеа Йои,, адареад ргойЕ, рогоиз /ауег, еачюпуагодупатис рагатее!, 
Беата сарасйу. 





* Тне гесеагсН {5 Чопе мил {Не #гате о! {Ие (пдерепаепЕ В&0. 
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